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Abstract 

The electronic structure of the doped manganites Pri-^Ca^MnOa was cal- 
culated by the LSDA+U method which takes into account the local Coulomb 
interaction between d-electrons of transition metal ions. In contrast to the 
standard local spin density approximation (LSDA) no Jahn- Teller lattice dis- 
tortions are needed to reproduce the experimentally observed orbital ordering 
in the undoped PrMnC>3. The correct ground state : a charge and orbital or- 
dered antiferromagnetic insulator, was obtained for Pr^Ca^MnOs. The 
results are in good agreement with the neutron diffraction data. 



I. INTRODUCTION 

The doped rare earth manganites Rei_ x A x Mn0 3 (Re is a rare earth such as La and A is a 
divalent element such as Sr or Ca) due to their peculiar correlation between magnetism and 
conductivity have been extensively studied during 1950s and 60s [1]]. The most throughly 
investigated was the Lai-^Sr^MnOs system. Undoped (x=0), LaMn03 is an antiferromag- 
netic insulator. Upon doping with Sr, this perovskite oxide becomes a ferromagnetic metal; 
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the connection between metallicity and ferromagnetism was well explained by the double 
exchange hopping mechanism ||. The discovery of colossal magnetoresistance phenomena 
in samples with Sr dopant densities in the 0.2< x <0.4 regime || brought a revival of the 
interest to these sytems. 

The Mn +3 ion in the hole-undoped compound LaMnOa has the high-spin d A electron con- 
figuration t\ ^e x g ^. The ti g orbitals hybridize with 02p orbitals much more weakly than the 
e g orbitals and can be regarded as forming the localized spin (5=3/2). In contrast to that, 
e g orbitals, which have lobes directed to the neighboring oxygen atoms, hybridize strongly 
with 02p producing in the result rather broad bands. The strong exchange interaction with 
t 2g i subshell leads to the splitting of the e g band into unoccupied e s j_ and half-occupied 
e g i subbands. The half filled e g ^ subband is a typical example of the Jahn- Teller system, 
and, indeed, LaMn03 has an orthorombic crystal structure [[| with distorted (elongated) 
MnOe octahedra. This cooperative Jahn- Teller effect is usually considered responsible for 
the opening of the gap in the half-filled e fl j band and an insulator ground state of LaMn03. 

The orthorombic crystal structure of LaMn03 (Pbnm space group) can be described as 
a perovskite with two types of distortion from a cubic structure: 1) Tilting (rotation) of 
the MnC>6 octahedra, so that Mn-O-Mn angles become less than 180°, and 2) Jahn- Teller 
distortion of the octahedra, with one long Mn-0 bond and two short bonds. With doping 
by Sr the orhtorombicity (the deviation of the values of the lattice parameters b and c/ y/2 
from a) becomes smaller and for x > 0.175 the stable crystal structure at room temperature 
becomes rombohedral, where only tilting of oxygen octahedra is present but all three Mn-0 
bonds are equal. Around the same value of the doping x, Lai_ x Sr x Mn0 3 becomes metallic. 

For Sr doping with 0.2< x <0.4, Lai-^Sr^MnOs is a ferromagnetic metal at low tem- 
peratures. However, as the temperature increases and approaches the Curie temperature T c 
a sudden increase of the resistivity is observed ||. As the temperature dependence of the 
resistivity p(T) has typically semiconducting behavior ( dp/dT < 0) above T c , this increase 
in resistivity at around T c is usually desribed as a metal-insulator transition. The colossal 
magnetoresistance effect happens when the temperature is close to T c : the external magnetic 
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field leads to suppression of the resistivity increase or to shifting the metal-insulator transi- 
tion to higher temperatures resulting in a very large negative value of \p(H) — p(0)]/p(0). 

The crucial point in the understanding of the effect of the colossal magnetoresistance 
is the nature of the metal-insulator transition with the temperature variation across T c . It 
was shown []5J that the double exchange mechanism alone is not enough to explain such 
transition. Millis et al. || suggested local Jahn- Teller distortions (Jahn- Teller polarons) 
as a main mechanism causing localization. Their idea is that a Ejt/W ratio (Ejt is self- 
trapping Jahn- Teller polaron energy and W is effective band width) is close to the critical 
value in Lai_ x Sr x Mn03 and as W decreases with destroying the ferromagnetic order above 
T c (the effective hopping takes the maximum value for a parallel spin alignment), this ratio 
becomes larger than the critical value to result in polaron localization. Varma [|7| argues 
that LaMn03 is not an insulator due to the Jahn- Teller distortion but is a Mott insulator 
and the Jahn- Teller distortion occurs parasitically. He explains localisation of the holes in 
the paramagnetic phase by spin polaron mechanism. 

By replacing La with other trivalent ions with smaller ionic radius, the Mn-O-Mn bond 
angle becomes smaller and the e g band width W is reduced to enhance the tendency to 
the carrier localization and lattice distortion. Variety of dramatic phenomena have been 
observed in (Pr, Nd, Sm)(Sr, Ca)Mn03 systems. As one of such examples, we study 
Pr!/ 2 Cai/2Mn03 in the present paper. This system has a very peculiar phase diagram 
||||. At low temperature, it is an antiferromagenetic insulator with a charge ordering of 
Mn 3+ and Mn 4+ accompanied by an orbital ordering: the orbital ordering and spin ordering 
is of CE-type. Between (~ 180 K) and T co (~ 240 K), there is no long-range magnetic 
ordering but the system remains to be insulating because of the persisitence of the charge 
ordering. Above T co , it is a paramagnetic insulator. 

The local spin density approximation (LSDA) |10[ was used by several groups for theo- 
retical investigation of the electronic structure of Lai_ x Sr x Mn03 |p"T|-|T3|1 . Those calculations 
have confirmed the importance of the Jahn- Teller distortion for correct description by the 
LSDA of the insulating antiferromagnetic ground state of the undoped LaMnC>3, because 

3 



the calculations with an undistorted cubic perovskite crystal structure produce a half-filled 
metallic e s j band. However in order to address the charge and orbital order observed in 
doped manganites it is necessary to go beyond LSDA by including the intra d-shell Coulomb 



interaction (LSDA+U method [p!4bfl6|). Recently such an approach was used for treating 
the charge ordering in magnetite Fe304, where LSDA gives a metallic uniform solution. 
The LSDA+U calculation gave a charge ordered insulator as the ground state of this com- 
pound [|17}]. In this paper we describe the results of the electronic structure calculation for 
Pri-zCa^MnOs for x = and 0.5 in the LSDA+U approximation. 

II. COMPUTATIONAL METHOD 

The problem of the charge ordering cannot be treated properly by the standard LSDA 
because of the non- vanishing self- interaction. In contrast to the Hartree-Fock approximation, 
where the self-interaction is explicitly excluded for every orbital, in the LSDA it is canceled 
to a good approximation in the total energy, but rather poorly in the local one-electron 
potentials, which are orbital-independent. The spurious self-interaction present in the LSDA 
leads to an increase in the Coulomb interaction when the distribution of the electron charge 
deviates from the uniform one. This effect can be illustrated in the following way. If we 
consider the system having less than one electron per site in the partially filled electronic 
subshell (for example e 9 -electrons in Lai-^Sr^MnOs) then the potential which this electron 
feels at some site must not depend on the orbital occupancy of this particular site, as the 
electron does not interact with itself. However, as the LSDA potential is a functional of the 
electron density only, increasing the electron density at one site and decreasing it at another 
one, with a formation of the charge ordering, will lead to an increase of the potential at 
the first site and a decrease at the second one. As a result, in the self-consistency loops, 
the charge distribution will return to the uniform density. In order to cure this deficiency, 



it is necessary to remove the self-interaction. The LSDA+U method p^JT6| is one of the 
practical methods to solve the self-interaction problem. 
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Let us introduce the d-orbitals occupation matrix defined by: 

1 f E F 



7T 



/ F ^ nlmMml (E)dE, (1) 



where G a , . , ,(E) = (inlma I (E — H) 1 I inlm a) are the elements of the Green function 
matrix (i: index of atom site, nl: main and orbital quantum numbers (for example 3d), m: 
index of a particular d-orbital and a: spin projection index). In terms of the elements of 
this occupation matrix {n a }, the LSDA+U functional [16|] is defined as follows: 



E LSDA+U [p a (r), {rf}] = £ L5DA [p CT (r)] + E u [{n a }} - E dc [{n rT }] (2) 

Where p a (r) is the charge density for spin-a electrons and E LSDA [p a (r)] is the standard 
LSDA (Local Spin-Density Approximation) functional, ^[{n^}] is the Hartree-Fock type 
electron-electron interaction energy given by 

E u [{n}} = |E{ m }, ff {(m,m" I V ee | m! , m'")n a mm ,n^J, ml „ 

(3) 

-((m,m" | V ee | to', to'"} - (m,m" \ V ee \ m'" , m'))n^ m ,n^„ m ,„}, 

where K e is the screened Coulomb interaction among the nl electrons. Finally, the last term 
in Eq.(0) corrects for double counting and is given by 

E dc {{n°}] = ~UN{N - 1) - ^ J[iV T (iV T - 1) + - \)\, (4) 

were N a = Trin^m') and N = + NK U and J are screened Coulomb and exchange 
parameters RI8| , |I9| . 

In addition to the usual LSDA potential, the variation of the functional (Q) gives an 
effective single-particle potential to be used in the single-particle Hamiltonian : 

H = H LS da + I inlma)V£ m , (inlm a | (5) 
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Knm' — Em"m'"{( m i m " I ^ee I m ' , m "') n rn/'m"' 

+ ((m,m" | V ee | m! ,m"') — (m, m" \ V ee | m'", m'))n^„ m ,„} (6) 

-U(N~l) + J{N° -\). 

The matrix elements of the screened Coulomb interaction V ee can be expressed in terms 
of complex spherical harmonics and effective Slater integrals F k GO 



as 



(m, m" | V ee | m', m'") = £ a k (m, m' , m", m"')F k 1 (7) 

k 

where < k < 21 and 

47T k 

a k (m,m',m",m"') = — — ^ (lm \ Y kq \ lm'){lm" \ Y k * q \ lw!") 

2fc + 1 q=-k 

For d-electrons one needs F°, F 2 and F A and these can be linked to the Coulomb- and 
Stoner parameters U and J obtained by the LSDA-supercell procedures [0,0 via U = F° 



and J = (F 2 + F 4 )/14 with the ratio F 2 /F 4 being to a good accuracy a constant ~ 0.625 



for the 3d elements ||T5|J21 



The LSDA+U approximation was realized in the framework of the linear muffin-tin 
orbital method in the atomic sphere approximation (LMTO-ASA) ||22|| . The values of [7=7.9 
eV and J=0.9 eV were used as obtained from the constrained LSDA supercell calculations. 

III. RESULTS AND DISCUSSION 

As a first step the electronic structure of the undoped PrMnOs was calculated. This 
compound has practically the same properties as LaMn03. The only difference is that due 
to the smaller ionic radius of Pr ion compared with that of La ion, the tilting of the oxygen 
octahedra in orthorombic Phnm crystal structure is stronger, resulting in smaller effective 
e g — e g hopping between Mn atoms and narrower band width. 
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The crystal structure of the doped Pri-^Ca^MnOa does not show transition to the rom- 
bohedral symmetry, as it is the case for Lax_a.Sra.MnO3, preserving the orthorombic Pbnm 
space group ||. However the values of the lattice parameters a, b and cj y/2 for x=0.5 are so 
close that it can be called a pseudocubic structure. In this structure only the tilting of the 
oxygen octahedra is present, but octahedra themselves are not Jahn- Teller distorted as in 
the undoped crystal structure. In order to study effects of the orbital polarization purely due 
to the intrashell d-d Coulomb interaction without any influence of the Jahn- Teller lattice 
distortion the calculation for the undoped PrMnOs was performed with the crystal structure 
parameters corresponding to Pri/2Cai/2Mn03. 



Kugel and Khomskii [£3] demonstrated for a wide range of Jahn- Teller magnetic com- 
pounds that the exchange mechanism (caused by an interplay between on-site Coulomb 
interaction U and intersite hopping t) by itself, without taking Janh- Teller distortions into 
consideration, is able to give a correct picture of the spin and orbital ordering. In this case 
the change in the lattice structure (the structural transition) is a secondary effect. One may 
say that the system shows the "Jahn- Teller ordering" without the "Jahn- Teller interaction". 
Actually an exchange mechanism acts jointly with the Jahn- Teller and the quadrupole mech- 
anisms. It is quite difficult to determine just which mechanism will be predominant in each 
specific case due to complexity of the real compound. One of the possible ways to answer 
this question in the frame of ab-initio band structure calculations is to use a crystal structure 
excluding distortions caused by the Jahn- Teller interaction. 



Such kind of approach was used for the cooperative Jahn- Teller system KCUF3 |L6 
It also has the perovskite crystal structure, but there is no tilting of the CuF 6 octahedra 
and all Cu-F-Cu angles are equal to 180°. Divalent copper ion in this compound has a d 9 
configuration {tlgi^lgi^^li) known to be a Jahn- Teller ion, and, indeed, the Jahn- Teller 
distortion is observed in the form of the elongated CuF 6 octahedra with long axes alterna- 



tively along a and b axes in the afe-plane. Kugel and Khomskii |^|] showed that a simple 
Hubbard model type calculation in a mean-field approximation without taking into account 
the electron-phonon interaction could predict the orbital ordering (alternating 3x 2 — r 2 and 



3y 2 — r 2 Cu3d-orbitals for electrons or x 2 — z 2 and y 2 — z 2 orbitals for holes) corresponding 
to the symmetry of the experimentally observed lattice distortions. 

The LSDA+U functional (eq.|2[) has a definite important improvement upon the LSDA: 
explicit cancellation of the self-interaction and the proper description of the lowering of 
Coulomb interaction energy with development of the orbital polarization. The functional 
(eqfj) is also rotationally invariant (the total energy does not change with the unitary 
transformation of the d-orbital basis set \m >), meaning that if the local coordinate system 
(and, correspodingly d-orbital basis set |m >) was chosen inappropriately for describing 
the resulting orbital polarisation, the unitary transformation which diagonalizes the non- 
diagonal occupation matrix rf mm , (eqfj]) will give a proper orbital basis set. 

The LSDA+U calculation for KCuF 3 in cubic perovskite crystal structure gave results 
fully confirming Kugel-Khomskii's prediction - alternating 3x 2 — r 2 and 3y 2 — r 2 orbital 
order. Moreover, it was not just the lowest total energy configuration among many possible 



stationary configurations as in the model calculation ||23|| , but the only one stable solution. 
The lattice relaxation to minimize the total energy is in good agreement with experimental 
data. Note that the standard LSDA did not show the lattice instability with respect to the 
Jahn- Teller distortion in KCUF3. We would like to emphasize that the orbital polarization 
caused by the d-d Coulomb interactions practically did not change with the lattice distortion 



24]. 

The use of the crystal structure parameters corresponding to Pr^Ca^MnOs (with 
undistorted oxygen octahedra) for PrMn03 may correspond to the use of the cubic perovskite 
crystal structure for KCUF3. However, for PrMnOs, we have to take into account the tilting 
of MnOe octahedra (present in Pr^Ca^MnOs structure) , which is essential for the proper 
description of the reduction in the effective e g — e g hopping due to the bending of the Mn-O- 
Mn bond. The LSDA+U calculation with the experimental A-type antiferromagnetic spin 
alignment resulted in a non-diagonal occupation matrix for spin density (for the e g subspace 
(3z 2 — r 2 and x 2 — y 2 orbitals) of one particular Mn atom): 



S 



0.31 0.22 



\ 



n 



1 



/ — n. 



(8) 



'mm' 



rn rn ' 



■ i — 



V 



0.22 0.58 



The diagonalization of this matrix gives two new e g ^ orbitals: 0i = 3y 2 — r 2 with occupancy 
0.71 and 02 = z 2 — x 2 with occupancy 0.18. For the second type of Mn atom 0i and 2 
orbitals can be obtained by transposition of x and y. 

The resulting orbital order can be presented graphically by plotting the angle distribution 
of the e 9 -electron spin density: 



As a Mn +3 ion has formally only one electron in partially filled e ff j subshell, the e s -spin 
density must correspond to the density of this electron (Fig.l). 

The plot of Fig.l and the form of the 0i and 2 orbitals diagonalizing the spin density 
occupation matrix show that the orbital order obtained in our LSDA+U calculations for 
PrMn0 3 is of the same symmetry as the one in KCuF 3 . The "parasitical" lattice distortion, 
which one could expect from such orbital ordering: elongation of Mn0 6 octahedra in the 
direction of the lobes of 0i orbitals, will reproduce the experimentally observed Jahn- Teller 
lattice distortion. 

The orbital polarization is far from being 100%, the difference in the occupancy of 0i and 
02 orbitals is only 0.53. The reason for this is the strong hybridization of e 9 |-orbitals with 
oxygen 2p-orbitals. In Fig. 2, the total and partial densities of states (DOS) for PrMn0 3 are 
shown (for e 9 -electrons DOS projected on X and 2 orbitals are also presented). 

The electronic structure obtained in our LSDA+U calculation is semiconducting with a 
band gap of 0.5 eV. Note that the standard LSDA needs the Jahn- Teller distorted crystal 
structure to reproduce the non-metallic ground state. The e^-band is split by this gap into 
two subbands: the occupied one with the predominantly 0i character, and the empty one 
with the 02 character. However there is significant admixture of 2 in the occupied subband 
and 0i in the empty one. 



p{0, 0) = £ {n, 



™L^ m (M)^(M) 



(9) 



mm 1 
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After checking that our calculation scheme gives a proper orbital order in the undoped 
PrMnC>3, we applyed it to Pr 1 / 2 Ca 1 / 2 Mn03. In order to treat the experimentally observed 
CE-type of antiferromagnetism it was necessary to quadruple the Pbnm-type unit cell and 
the supercell had 16 formula units. The distribution of the Pr and Ca atoms in the lattice 
was chosen according to a model where all Pr ions have only Ca as their nearst neighbores in 
the lanthanide sublattice. We did not impose any symmetry restrictions in our calculation, 
the integration in the reciprocal space being performed for the whole Brillouin zone, and the 
self-consistency iteration was started from a uniform distribution of e 9 -electrons over both 
e 9 -orbitals of all Mn atoms. 

The result of the self-consistent calculation was a charge and orbital ordered insulator. 
The e 9 -electron spin-density-plot calculated by eq.|9| is presented in Fig. 3. This is in a striking 
agreement with the orbital order derived in || from the neutron diffraction mesuarements 
(Fig. 4). The unexpected result is that the total number of d-electrons at all types of Mn 
sites are nearly equal (4.99 and 5.01), so that formally Mn 4+ ions (in the Fig.3 the ones with 
symmetric in-plane density distribution) and Mn 3+ ions (the ones with the density strongly 
anisotropic) have nearly the same number of 3d-electrons. However the difference in the 
magnetic moment values is more pronounced : 3.34/i# for Mn 4+ and 3.44 Hb for Mn 3+ . 

IV. CONCLUSION 

The results of our LSDA+U calculations for Pri-^Ca^MnOs demonstrate that taking into 
account the Coulomb interaction inside the d-shell of Mn together with the intersite hopping 
is enough to reproduce not only the orbital order present in the undoped manganites, but 
also the localization and ordering of the holes in the doped materials. Naturally the lattice 
relaxation (Jahn- Teller and polaron lattice distortions) will significantly renormalize the 
total energies, gap values, charge disproportion and magnetic moments. However it seems 
that the symmetry of the ordered state can be obtained without explicitly including the 
electron-lattice interactions. This is in clear contrast to the standard LSDA where energy 
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band splitting and orbital polarization appear only as results of the Jahn- Teller crystal 
structure distortion. 
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Figure captions 

Fig. 1. The calculated angle distribution of the e^-electron spin density in PrMnOs. 
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Fig. 2. The total and partial densities of states (DOS) from the result of LSDA+U 
calculation of PrMnC>3. a)Total DOS per unit cell (4 formula units); b)Partial e g DOS of 
Mn. Solid lines are for the 0i-projected DOS, and dashed lines for the 02-projected DOS 
with <pi and </>2 -orbitals diagonalizing the occupation matrix.; c) and d) Partial DOS for 
two types of oxygen orbitals. 

Fig. 3. The calculated angle distribution of the e 9 |-electron spin density in 
Pr 1/2 Ca 1/2 Mn0 3 . 

Fig. 4. The scheme of the spin, charge and orbital order in Pri/2Cai/2Mn03 deduced from 
the neutron difraction data ||. The open circles with the lobe of the e g electron density 
distribution denote Mn 3+ and the filled circles Mn 4+ . The arrows denote the magnetic 
moments. 
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